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ABSTRACT 


The nai& objectives of this study is to develop s computer code 
that can assess the performance of the thermal parameters of the 
first wall during the steady state normal operation of the 
"ADITYA" (The Indian tokamak). 

A thermal analysis model is presented on the basis of which a 
code is developed. The code is capable of calculating the j 

I 

temperature profile in the first wall. Host of the parameters' 

1 

required for these calculation are inbuilt within the code. Two | 

types of density and temperature profile of electron and ion are < 

1 

considered. 

Poloidal distribution of the thermal flux (Due to , 
bremsstrahlung, cyclotron, charge exchange, line radiation, and 
neutron) are calculated, with the assumption of a specific coolant 

I 

arrangement system this is then coupled to a one dimensional heat- i 
conduction equation to derive the temperature profile in the first i 

i 

wall. I 
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CHAPTER 1 

IMTRODUCTION CAND HRST WALL PERFORMANCE) 


1.1 Introductrlom 

Th«rmon«cl«*r fusion rosoarch is nsklnt rapid proftrsss 
around ths world. Tokaaaks havs so far ylsldsd bast rosults 
anongst ths various nafinstie eonflnsnant sehaass undar study. Tha 
Joint Europaan Torus (JET), has racantly raportad conditions vary 
closa to thosa raquirad for braak-avan. Pralinlnary dasign 
stttdias hava alraady baan eoatplatad for tha naxt ganaration 
tokanak eallad ITEE (Intarnational Tharaonuclaar Exparinantal 
Raactor). Fusion raaetors siaant for alactric powar tand to ba big 
and axpansiva. Tka first coaaarclal application of tokasudcs in 
India is likaly to ba as a low gain sourca of fusion-nautrons for 
flsslla fual production. Study of a fusion braadar is underway 
at tha Institute for Plasaa Rasaarch (ZPR), Bhat, Gandhinagar 
(India). [1] Objactivas included in tha progranna are, 

datarninatlon of tha tokanak and blanket paramatars raquirad for 
tha device and a study of how these devices could fit into tha 
Indian nuclear powar prograsuae. Tha prograaoia at IPR. is divided 
into two types, l.a. short-tarn and long tarn objactivas. In tha 
short tarn ob jaetiva nain - anphasls is on a low fissile - 
production rata (FPR « few kg/yaar) and reducing cost, 

design. Tha device is intended to danonstrata tha 
feasibility of fusion breeders. Tha long tarn objaetiva is a 
*cala-up of tha first. It will hava a nodarata FPR (few- hundred 
kg/yaar). Its larger, sisa and nora eoaq>lax design will incraasa 
tha energy gain to a level that is anargatically attractive. 
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Soae l&tttrttedl&tft steps sre necessary between tbe currently 
existing Indian tokanak *AD1T7A' and the final device. llork Is , 

In progress for deternlnlng the opt Inal parameters of the 
Intermediate steps. 

1.2 Tdkamak and First Wall 
Tokamak [2-5] 

If a nucleus of a deuterium fuses with a nucleus of a 
tritium an «t-partlcle Is produced and a neutron Is released. The 
nuclear rearrangements results In a reduction of total mass and a 

I 

consequent release of energy In the form of kinetic energy of the i 
reaction products. 

! 

2^41 
•D + jT > jHe + ^n 

(3.5 HeV) (14.5 HeV) 

In order to produce the fusion of nuclei of deuterium and ! 
tritium It Is necessary to overcome the mutual repulsion due to | 
their positive charges. Fusion brought about by heating the ; 

j 

nuclei to a high temperature Is called thermonuclear reaction. 

i 

t 

The temperature required to start fusion Is approximately of the | 

I 

A ' 

order of 10 KeV (•» 10 K). At such temperature the fuel (Da T) 
fully Ionised. The neutral gas thus formed Is called plasma. 

I 

A tokamak Is a device, where this plasma particles are | 

t 

f 

confined to a toroidal region by a magnetic field. In a tokamak | 
It la necessary to confine the energy of a sufficiently dense 
plasma for a time which allows an adequate fraction of the fuel | 
to react. | 

f 

In a tokaxMk at the required ion-density of around 10^** 
this eohflnement time should be around 1.5 sec. [5]. The product 
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o£ ion density and confineaent time is often knovm as Lawson 

I 

I 

criteria, which says, the energy output froa the fusion will be 

greater than energy input if the product of ion-density and i 

20 -3 

confineaent tiae is around (1.5-3)xl0 a sec. (vthen teaperature , 
is around 10 20 XeV) [4]. 

The principal aagnetic field in a tokaaak is the toroidal 
field, however this field alone does not allow confineaent of , 

plasaa. In order to have an equilibriua between the plasM 

pressure and aagnetic pressure, it is necessary to provide a ! 

poloidal field. Poloidal field coils are used for this purpose, | 

but aainly the plasaa current contributes to the poloidal field. | 

I 

The layout of principal coaponents in a conceptual tokasAk > 
reactor is shown in Fig. 1.1 and 1.2. 

i 

The core is aagnetlcally confined plasaa. Moving radially | 

i 

outward the next feature is the liaiter or a divertor (not shown 
in Fig). The aain use of liaiter is to scrape off and reaove 
impurities and fuel ions froa the surface of the plasaa as they i 

i 

aove outward but before they can hit and damage the first wall. ' 

Diverter plays the same role as the liaiter, the difference is j 

that diverters are aagnetic liaiter, there exist a diverter ' 

I 

current to produce suitable aagnetic field which guides the { 

, ' i 

escaping plassMi particles to a - separate chamber where the j 
particles are neutralized and pumped away. An annular insulating 
vacuum region separates the plasaa froa the first wall. The 
first wall(soaetiae knoim as the vacuum wall) is the material 
surface that faces the plasaa and through which the neutrons and 
also the radiant energy flux passes. Behind the first wall are 
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coolant paaaafioa daaisnad to raaiova tha thornal boat loading £roM| 
radiation# and particle bonbardnant. Bayond tha firat wall ia a 
blanket, which in D-T raactora contain# lithiun to bread the; 
required tritiun fuel. Ih fiaaile fuel breeder type tokaaiak, 
thia blanket contain# auitable auiterial, which after reactinft 
with the fuaion neutron# produce# fiaaile naterial. 

jLi + Jn ► Jse + ♦ 4.8 HeV 

»L 1 . , Jh. ♦ It . Jn" - 2.5 H.V | 

Blanket thickneaaea of around i a. are required to noderate and | 


abaorb 14 HeV neutron# and to produce required amount of tritiun 


or fiaaile material. Typical breedinft reaction# are: 


238 

92 


U 


+ 


239 

94 


Pu 


232,.. ^ 1_ 2 /9 . 233 

^jjTh + ^n ► ,2 


U 


Beyond the blanket ia a nasnetic thermal ahield, which alow# 
down and abaorba neutron#. The auperconductinq aasnetic coil# 
are next major ayateaa in the reactor. Outaide the 
auperconductins ma£netic coil# a biological ahield ia preaent to 
reduce radiation to low enough level that maintenance peraonnel 
can approach the exterior of the reactor. 


First Wall [2>&} 

firat wall of a tokamak reactor and it# aaaociated 
atructurea operate in a aevere environment like 

1. Neutron- r el eaaed in the fuaion reaction and the charge 
exchanged neutral# the wall 
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2. Th« vail la aubjactad to a vary hifih praaaura 
diffaranca» i.a. tha aida facinc plaaaa la at vacttuai and othar 
a Ida at nornal praaaura. 

3. Haatlnft fron ona and (aida faeln^t tha plaaaMi) and 
coolii^ of tha outar aida. 

Baeauaa tha Ufa of flrat wall dapanda on fluanca of 
nautrona, tharaal anargy and charsad particlea In a tokaaak, and 
thara ia a atronfi notlvatlon to incraaaa thaaa fluxaa to their 
naxinun poaaibla valuaa, ao thara ia a naceaaity of raachlng an ; 
optinun balance batvaen allowable wall loadinsa and firat wall i 
llfatiaa. 

Tha firat wall conaiata of a chanbar that containa tha j 
plaana and tha aaaociatad atrueturaa. The aaaociatad atructuraa i 
include linitara or divertara, and low Z-natarial coatinfi. Tha 1 
dlvertara parturba tha confinins nasnatic field* in a way ao that ' 
tha aacapin& fuel iona inpinga in a controlled way on tha 

1 

diverter plataa* uaually remote from tha firat wall. A vide 
variety of firat wall material haa bean reported in literature. 
Tha moatly uaad material ia sanarally S.S. *- 316 or S.S. - 314 

'type. Tha exact dimanalona of firat wall dapanda on othar 
tokamak paramatara* however tha thicknaaa llaa batvaen 1 cm-Z cm. 

Tha functiona of tha firat wall are; 

!• To remove heat cauaad by bramaatrahluns and aynchroton or 
cyclotron radiation, by fuaion nautrona and by fuaion 
reaction producta which bombarda tha firat wall. 

It copaa with tha charged partic-laa and charge— exchanged 
anargatic particle fluxaa that ariaa . in tha. plasma and 


2 . 
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deposited on th* first vsll. 

3. Ths first wall dsals with loeslizsd snsr^y dspositlon 
rssaltinfi fro* plasaa disruptions or run away eloctrons. 
Bscauss of thsss proesasss ths plasaa duaq;>s a sianif leant 
fraction of ths plasna snsrey on ths surface of ths first 
wall . 

i. Ths first wall designed in a aanner so as to reflect back 

most of ths cyclotron radiation into ths plstsma, hence 
prevent cooling of plasna rapidly. 

1.3 system Cods Tor Tokamak 171 

A computer code is needed for integrated modelling of the 
tokamak. Such a system code is required for the following 
reasons : 

1. A tokamak based breeder has a large number of free 
parameters, along with a number of constraints. It is necessary 
to determine a wide range of optimal parameter sets from the 
point of view of cost as well as required parameters of 
sub-system. This can only be done by a system code. 

2. At different stages of the design it may be necessary to 
introduce new information and constraints. The code will then be 
required to determine the revised optimal parameters. 

A systems code allows a wide range of parameter sets to be 
*3camined within a reasonable time frame, and the optiauil set to 
be selected. Estimates of system costs are^ however only 
approximate. On the other hand, a detailed design gives accurate 
values for the cost and parameter values, but is . expensive in 
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t«ra8 of coa«>utotional ti*« and offort. Henco two approaches are 
coapl enentary . 

1.4 Earlier Work 

Systens codes already exist, and have been used extensively 
for conceptual desifins of aany tokaaaks through out the world. 
Several laboratories have indicated that they are unable to 
transfer their codes due to variety of reasons [7]. Hence the 
need for developing a code. 

A nulti code consisting essentially of three nodules, the 
'first wall*, the 'neutral gas* and the 'plasaa* stodule is 
developed by A. Nicolai and D. Reiter [8]. The first wall nodule 
contains the thernal aspects of the first wall; on the basis of 
inplantation of charge exchange neutrals and ohnic heating only. 

K.C. Carroll and 6.H. Hiley have reported a nodel for first 
wall thernal response [9}. They have not considered the effect 
of line radiation and radiative reconbination radiations fron the 
inpurlties present in the plasna as a source of heating of first 
wall. 

H.C. Carroll in his Ph.D. theses [10} developed a very 
accurate nethod to conpute the local wall power loading, using 
the toroidal correction factors. A polynonial fit fornula with 
data are available to correct the local wall power loading for 
toroidal curvature. 

1.5 ^e»snnt Ncur-k 

The naln objective of this study is to develop a conputer 
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cod« that can asaeas tha parforaance of tharaal paraaatara of tha 
fXrat vail during noraal oparatlon of tha ’XDITYA* tohaaak Tha 
tara 'noraal* danotaa conditlona undar which plaaaa paraaatara 
vary in a controllad aannar. Thia ia diffarant froa tha plaaaa 
'diaruption* whara tha anar^y of tha plaaaa la duapad at tha 
firat wall in a vary abort parlod of tiaa. Noraal oparatlon ia 
hovavar not raatrictad to a ataady-atata oparatlon, but tha coda 
davalopad hara ia aaattainfi a ataady-atata noraal oparatlon only. 

A tharaal analyaia nodal ia praaantad on tha baaia of which 
a coda ia davalopad. tilth aoaa fundaaantal tokanaJc paraaatara aa 
input, tha coda can giva tha taaparatura in tha firat wall - aa a 
function of (x, 0) , whara x ia tha dapth into tha firat wall 
natarial and *d* ia tha poloidal ancla. Noat of tha paraaatara, 
lika varioua raaction rataa, ara inbuilt within tha coda, and for 
diffarant input conditlona, thaaa raaction rataa ara 
autoaatically avaluatad froa aui table polynonial fita diacuaaad 
latar. 

Tha coda can ba uaad for two typaa of profilaa of danaity 
and taaparatura, i.a. flat type and paafcad type. 

Plaaaa conaidarad inaida tha tokanak ia of Dautariun and 
Tritiun 50t nixtura. It ia vary difficult to obtain a 100% pure 
aixtura of D & T for thia purpoaa. Hydro^an ia a coamon iapurity 
in both DAT. Alao bacauaa of aroaion froa tha firat wall which 
ia autda of S.S. - 314 (in aoaa caaaa other ^rada of S.S. ara 
uaad) which haa high carbon contanta aatarial, ai£nif leant anount 
of carbon iapuritiaa alao liaa within the plaaaa. It haa bean 
conaidarad hara that carbon and hydro£an ara tha two iapuritiaa 
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in the plasaa of D & T. Existence of other inpurity atoms 
however can not be ruled out, but to make the model simpler, only 
two types of impurities are considered. Since the code developed 
is in a modular fashion, inclusion of other impurities if 
necessary in future can be done with some minor modifications. 

The torus cross section considered here is elliptical one, 
which is the case of "ADITTA”. First wall material chosen is 
S.S. - 314. A specific combination of coolinfi arransement of 
first wall is chosen i.e. water flowinfi in tubes attached to the 
back of the first wall, with one/two input and one/two output is 
considered. 

This code can however be used for other type of torus 
cross-sectional seometry and coolant arranfiement, because of its 


modular fashion. 



CHAPTER 2 

THEORETICAL MODEL 


2.1 A Brief Outline of the Method Followed 

In this chapter nethoda and aquationa are developed to 

execute the calculation aequence diacuaaed in 2.3 and 2.4. 

Theoretical development not diacuaaed in other referencea ia 
deacribed in more detail. Uhere other referencea [3-5] are 

involved, material la briefly reviewed for the aake of 
completeneaa. To adequately explain aome of the numerical 
procedurea, it vaa neceaaary to include calculational examplea, 
which are included for "ADITTA” parametera uains thia model in 
Chapter 4. 

A thermal analyaia model, with correapondins aolution 
methoda haa been developed for determinins temperaturea at the 
aurface and interior of the firat wall. Given the wall material 
and £eometry, thia model ia capable of computins the local wall 
power loadinga (including cyclotron, bremaatrahluns, charge 
exchange, line and radiative recomb Inationa radiatlona, and 
fuaion neutrona). The model encompaaaea a geometric analyaia for 
determining loading variatlona along the pololdal anglea 
coupled with a one dimenaional heat tranafer analyaia for 

determining temperaturea in the firat wall. 

The two typea of profile conaidered are flat type and peaked 
type. In a flat type of profile, the denaitiea and temperature 
remaina conatant in the entire plaama croaa aection in aide the 


torua 
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C2.1) 

C2.2) 


n, T, is constant. 

* » s * » * 

The psaked type profile is asstmed as: 

where a. « major radius of the plasma column. 
oL^r are constant. 

The heat transfer problem analyzed is pictured in Fift. 2.1 
[10-11]. Radiations, charge exchanged neutrals and fusion 

neutrons are emitted from the plasma. After passing through a 
vacuum, they impinge on the first wall. Charge exchanged 
neutrals, cyclotron radiation and line, radiative recombinations 
radiations deposit their energy at the surface of the wall, 
whereas bremsstrahlung and neutrons penetrate into the material 
of first wall; in the case of neutrons, it may pass through 
completely the wall. Heat from the energy deposition is 
conducted to the bulk of the coolant and is then carried out of 
the reactor to be converted to useful energy. So, the source 
terms, which are responsible for heating of first wall, is 
divided into two types, i.e. surface heating terms and volumetric 
^^••ting terms. They are then treated separately for the 
temperature solution. 


The sequence followed in this model is shown in block 


diagram in Fig. 2.2. 
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— »A&XATXOM~> 

r/m. 

Bremssirahlung 


***r*»**jjr*r»j*** •!s 



Fifi 2.1 


Sketch of the FIrst-v/aU Heat 
Transfer Problem. 




Fig 2.2 

Calculation Sequence 




2.2 DvtAils of Soiuro* Torae 

Different proces 4 i «8 in th« plaaaut which contribute 
sifinif icantly to the well loedin£ are 

1. Brensetrahlunfi radiation 

2. Cyclotron radiation 

3. Neutrals foraation due to charge exchange nechanisn 

4. Line radiations 

5. Fusion neutrons. 

Details of these source terns are discussed below, with equations 
and procedure for their calculations. 

2. 2. 1 Brensstrahlung Radiation 

This forn of radiations occurs in a plasna when one 
particle encounters other and aake a sudden change in direction. 
It is generally identified as a free-free transition radiation 
[12]. It constitutes an energy loss and cooling nechanism for 
the plasna which increases with the increase of proton nunber of 
the fuel species or inpurities. Brensstrahlung radiation is 
enitted in a continuous distribution of frequencies ranging fron 
radio frequencies upto photons with energies several tines the 
electron kinetic tenperature. The average brensstrahlung photon 
cones out of a plasna with an energy conparable to electron 
tenperature of the plasna [3]. 

Since the electron tenperature of a fusion plasna renains 
sufficiently high (> 10 KeV - 10*K) there is a negligible 

probability that these brensstrahlung photons will be reabsorbed 
within the plasna. Such energetic photons are - absorbed by the 
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first wall without sifinificant reflsction from it. 

A ri£orott8 formulation for the bremsstrahluns radiation 
power loss can be developed by using quantum mechanics and 
treatment of relativistic effects [13]. But here, an approximate 
analysis and in a non-relativist ic range is considered for the 
formulation of bremsstrahlung radiation. 

From [3], total power 'p* in watts radiated from a single 
electron of charge 'e* moving with a velocity V and acceleration 

m 

V is given by the formula: 


dF 

WE 


2 


6n « c' 
o 


[l Y. - <V X V)^/c^ 
^ Cl - V^/c^)^ 


1 „ 


(2.3) 


In the non relativistic case 


WE - 


e 



6n 


U 


(2.4) 


Assuming an electron of charge 'e* is approaching a positive 
ion of charge 'Ze' with impact parameter 'b* and during the 
interaction the electron suffers a force which results in an 
emission of a photon of energy h 2 > (Fig. 2.3). 

Average Coulomb force between the electron and the ion when 
separated by a distance '*b* impact parameter is 


*e • ^ 


Ze' 


4n « b 
o 


(2.5) 


Average acceleration of the electron, when it interacts with 
a ion of proton ho. 'Z* is from (2.5) 
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Electron gyrattnr in a atronc nafnetie inducUoa and amlttinf ayachrotroa 
n^tloa alone its rolocity roctor la ttw plana of Its orbit. 


Fifi 2.4 
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Z« 

I 


2 


(2.«) 


4it b 

O e 

Now (2.3) can b« aiodlfiad as. 


dp _ 

^ " ox 3 3 3 J2 ri" 

96 R « c 1 B_ b 
o e 


C2.7) 


If the electron was in the vicinity of the ion for a tine 't', 

t = b/V (2.8) 

Radiation per collision is 


B « t X ^ (2.9) 

««> ^ ^ Joules (2.10) 

In unit tine within inpact - paraneter (b) to (b+db) the 
nunber of collisions in a unit volune between electrons and ions 

t sec = n^ X n^ X V X 2fib db (2.11) 

Total power enitted by the electrons in that tine is fron 

( 2 . 10 ) 

dp = 2nb^ nj n^ || db U/cn^ (2.12) 


Substitutine ^ fron (2.7) 


dp 


2nb n.n 
1 e 


« 

O B 


db N/cn' 


(2.13) 


Total brensstrahluns power 
obtained by intesratinft (2.13), 
to b i.e. Oft. 


enitted in a. unit volune can be 
fron b^i^ i.e. conpton wavelenftth 


The conpton wavelencth is siven by 
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'nin 


h/2« 


Total power emitted from a unit volume of plasma is 


Z^e^ n, a V 

P = - y * W/cm 

24r o'* m h 

o ^ 


(2.14) 


Substitutin£ the characteristic electron 
thermal velocity. 


V 


8 K 1/2 

« ** J 


velocity 


'V* by mean 


(2.15) 


P = 


»2 6 

2 e n. n 
X e 

24r c^m h 

O 0 


8 K T 


1/2 




W/cm* 


(2.16) 


which sives P = 1.69xl0”^^ n^ n^. W/cm^ (2.17) 

where is in eV, densities are in # cm"^. 

The total bremsstrahlune radiation from the plasma when 
there are ions of different species can be eiven by. 


P = 1.69x10"^^ n^ 2 “i 2^ W/cm^ (2.18) 

i 

Calculation of ^ **i^i * 

1 

In the plasma, carbon (impurity) exist in 7 different 
states. Hydrofien (impurity). DAT (fuel) exist in 2 different 
states. 

In Ref. 14 Post et al have reported polynomial fit- data for 

47 elements to compute S n, zf, as a function of electron 

i - 

temperature (0,002 KeV S T^ S 100 KeV). These data are senerated 
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U8in£ a "coronal aquilibrittn* nodal* i.a. tha disbribubion of 
char«a atata of any ion is Aovarnad by tha aquilibriun betwaen 
alactron inpact ioniaation amd raconbi nation. 

Tha polynonial fit is sivan by: 


- 5 r ,i 

< > = 2 n. * 2^ c(i) [lofijj, T^(KaV)J (2.19) 

cCi) ara availabla for 47 alanants (26 £ Z ^ 92) in [14} in a 
tanparatura ranfia of (0.002 KaV £ £ 100 KaV). 

\ nj^ for carbon is calculated usins this fit. 

To calculate y ^i^i hydroftan, a£ain coronal 

aquilibriun nodal is considered. Electron inpact ionisation 
(Eli) for hydro£arn is 

a^ + H(1S) + E* + 9^ 

and tha radiative raconbination (RR) for hydro£an is as follows: 

a + H* 1 - H(1S) + hw. 

Tha rata of these reactions ara 

< O' • V >£jj ^(IS) < O' . V n^ n ^ for RR. 

H 

A88unin£ coronal aquilibriun* i.a. an aquilibriun between Eli and 
RR [16] 


< O' . V > 


Eli % *^(1S) 


< » . y >„ n. 


( 2 . 20 ) 


< . V > 


RR 


< ^ . V >g„ E 


( 2 . 21 ) 


“H(1S) 

Further since n^ » total hydro£an density is the sun of 
ng(lS) and n ^ or 


•h 


“H(1S) * 


( 2 . 22 ) 
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“h 


fi ♦ 1 

I < . V J 


Eli 


and 


“H(1S) “ “ *^g+^ 


(2.23) 


(2.24) 


< a . V >£j 2 fox* hydrosen is calculated fron polynoaial fits in 


reference [15] 


In < o- . V > = 


= 2^ 1>„ in CT,)“ 


(2.25) 


b are available in [15]» for a temperature ranse 0.1 eV to 20 
n 


V for hydrossn is calculated fron the 


KeV. < a 
seai-eapirical relation in [15]. 
For n s nl s IS; 2S; 2P; 

< . V 


nx xs; xo; xr; . 

>R* - »nl * ‘® L®a '*^\ ~Tn 

* e 


cm^/sec 


n nl 


(2.26) 


where Ry 

s 

Rydbers 

constant » 

13.58 eV 


T 

= 

Haxwell 

temperature. 



^n 

s 

/ T 






nl 

IS 

2S 

2P 



^nl 

3.92 

2.47 

6.22 



*nl 

0.35 

0.12 

0.61 

For n 2 3 






< «- .V > 

£R 

= 5.201 X 

10-“ 

U 

®1 ^n^ 

3 

(/9^) cm /sec. 


(2.27) 
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E^C/3^) B exponent ial intesrel. 

In the code developed only transition from 2S to IS is 
considered. 

So now, (2.23) and (2.24) is solved to obtain population of 

2 

different staftes of hydro£en froa where £ n.z, for hydrogen is 

i ^ ^ 

obtained . 

2 

Calculation of S n. s. for D & T is done assuninfi no 
i ^ ^ 

existence of D & T neutrals in the plasna. This is because the 

14 

total number of D & T is very high of the order of 10 (for 
'ADITYA* design), in conparision to to D^, (neutrals) at 

high tenperatures (> 10 KeV) [14.]. So to save computation time 
without any significant error, in the calculation here existence 
of D & T neutrals is neglected in the code. 

C 3 rcloitron Radiation 13,171 

Uhile bremsstrahlung is a process that takes place within 
the short time an electron spends in the vicinity of an ion or 
atom, cyclotron emission is a non-collisional mechanism that 
occurs in the interval - between collisions. It originates from 
the orbital acceleration of an electron in a magnetic field. 
Cyclotron emission by ions is insignificant because of their 
large mass . 

For electrons whose energies are less than several hundred 
electron volts, cyclotron emission is concentrated in a single 
spectral line, whose radiation frequency a equals the orbital 
frequency, (u^^ b eB/m, B » strength of magnetic field). 
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A £yratinfi electron in a etrons nasnetic field (Fifi.Z.d) 
with characteristic frequency '<»* emits photon because of 
centripetal acceleration [13]. The radiated power of a 
relativistic electron is fiiven by [3],, equation (2.3) 


<i» 



(1 - V^/c^)^^^ 


(2.28) 


In terns of this relativistic syrofrequency, the velocity and 

acceleration of the electron may be written as 

V « «V (2.29) 

■ • 
where V- acceleration, for fiyro motion is perpendicular to V. So 

for a sinsle electron from (2.3), 


dp 

Ht 


2 „2 

e u 


<n c"' (1 - y^/c^) 


watts 


(2.30) 


Substituting for a in (2.30) 


dp 

3t 


e*vi 


watts 


(2.31) 


Total power radiated by all electrons from unit volume of plasma 
is, n times the radiated power in eq. (2.31). The electron in a 
plasma of fusion interest will have a Haxwellian distribution; of 
perpendicular velocities. Integratinfi over such a Kaxwellian 
distribution from [3]. 



ip 


e" b 2 n. 


K T 


3n « 


(m^c 



5K T 

e 


. .J W/cm^ 
(2.32) 



(2.33) 


This can be approximated to [18]. 

P « 6.21 X lO"*® B® n T U/c«® 

^ B 

where B la In Caues 

n^la in # cm"® 
in eV. 

The photons which are emitted in cyclotron radiation from 

fusion plasma have very high frequencies approximately 100 GHs 

14 -3 

[3]. In fusion condition plasma densities (10 cm ) and 
temperatures (10 o' 20 KeV) such high frequencies photons are 
reabsorbed. Such photons are also perfectly reflected from the 
first wall. Reflectivities at such high frequencies are nearly 
0.98 or even higher [3]. Since most of the cyclotron radiations 
are absorbed in the plasma and because of a high reflectivity, it 
is very difficult to estimate the exact amount that escapes to 
the wall. Here in the present model it is assumed that only 2 
percent of the total cyclotron radiation is reaching the wall. 

2. 2. 3 Charge Exchange MeoLrals 

Charge exchange is a process where redistribution of 
particles takes place 

A + B* ► A* + B 

Here an electron goes from atom'A* to the ion 'B^*. The 
electron transition from the atom "A* to the ion 'B^* is 
accompanied by transfer of the momentum *mv*, where m - electron 
mass and v s relative velocity of the nuclei of 'A* and 'B**^*. 
This process occurs in a plasaut and leads to significant loss of 
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energy from the pl&ama. The charge exchange neutrals nove quite 
fast and dump their energy at the vail surface. 

Out of the various charge exchange reactions possible in the 
plasma, some significant contributions comes from the reaction 

—+ m.O ^o ..+ 

H + H *■ H + H 

The rate of other reactions are very small [19]. So only the 
above reaction is considered as a source of charge exchange 
neutrals formation. 

Kate of the above reaction is given as 

»cx - V ^ 

The factor n ^ and n calculations are described in detail in 
H* 

section 2.2.1. 

From [IS], polynomial fits are available to compute 
as a function of temperature T of the plasma particle 
and energy E of the incident particle 

8 8 

ln<«'.V> “ 2 2 *^lj (2.35) 

i»0 j«0 

The coefficient a^j are available in [15]. These fits are 
valid in the range 0.1 eV < E < 20 KeV and 0.1 eV < T < 20 XeV. 
Energy of these charge exchanged neutrals are very high 
approximately 1.5 times the ion-temperature in the plasma. 

So in terms of energy, the amount generated from unit 
volume of plasma, because of CZ neutrals is 

^CX * “o V < . (1.5) T^ H/cm^ (2.36) 

where T^^ = 'Joules. 
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2. 2. 4 LIm Radiat.lonB 

Atonic line radiations can be divided into three 
categories [3]: 


(i) Ordinary spectral line radiation 

Ordinary spectral line radiation occurs when an 
electron is excited to a higher state by electron impact. It 

-7 

usually de-excites in a very short time 10 sec by the process 
of an emission of a photon. 

* 

,++ ,++ ,++ 

A + e ► A + e ► A + hi* + e 

(11) Radiative recombination : 



Radiative recombination occurs 

as a 

result 

of 

a 

capture of 

a free electron by an 

ion with 

a third 

body 

in 

the 

vicinity. 

Binding energy of the 

electron 

is radiated 

when 

it 


falls to ground state 

* 

, «++ » «+ 
A + B + e 1 A + B 

* 


B^ 


-♦ B + hi* 


(iii) Dielectronic Recombination 

Dielectronic recombination involves capture 
electron with excitation of a orbital electron at the 


such that both electrons are in excited state 


«++ 

A + e 


A 


*** 


U 


-♦A + hi»w 


+ + 

■* h e 


There are two ways in which the doubly excited 
decay. 


of a free 
same time, 

+ 

atom can 


(i) either it will give rise to two spectral lines 


(dielectronic recombination) or 
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(il) ionized with the release of an electron 

(autoionization) 

The autoionization process is very fast on the order of 

— 13 

10 St so only dielectronic reconbination is of interest. 

Atomic line radiation can not be emitted from completely 
stripped ions of the light element fuels of fusion Interest; but 
it can occur from incomplete ionized or neutral of high Z 
impurity. 

In this model, line radiations from impurities (carbon and 
hydrogen) is only considered. 

Line radiations from hydrogen: 

In case of hydrogen, out of the three processes discussed 
above, spectral line emission and radiative recombination is 
possible. 

Spectral line emission: 

Reference [18] gives empirical relation for spectral line 
emission radiation estimation. 

For a transition from n *■ m. 


= 5.1 X 10 

ziin 


-25 ^nm *1 ^^e 


g T 
•m e 


n n j. AE ^3 ^ AE 


U/cm” 
(2.37) 


where f^„ = oscillator strength, 

ZlHI 


th 


g^ - statistical weight of m state (m=l, ground state) 
T in eV 


n^ = neutral density 

AE = E - E 

nm n m 


To make the model simpler only transitions to ground state (m^l) 
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&re considered front n ^ 2 to 6 only. 

So 

T - R i in-25 ^nl % % 

^nl " 5.1x10 -JJ 2 exp 




(2.38) 


The oscillator strength (n = 2 to 6) are available In 


[15] and frost [18], 


= 13.6 




(2.39) 


Radiative recostblnatlon: 

The radiative recostblnatlon process Is as follows: 

.+ 


e H 


E(n) + hu 


n :£ <e 


only ground state (n=l) Is considered In this stodel. 

Rate of this reaction Is 

R = n n^ <o'.V> 

® H”" 


(2.40) 


Calculation of <o’.V> for this reaction Is given In (2.26). 
Further energy of the photon released In this reaction Is 

Ry 


E. 


hi» 


^ ^ 
n 


(2.41) 


Ry = Rydberg constant = 13.58 eV. 

Energy released from unit voluste of plassta through 
process is thus. 


E = n n . <o'.V> 




+ ^ 1 W/csi^ 

Tk^ J 


this 


(2.42) 


Line radiations from carbon: 

Post et al in Ref. [14] have given polynostlal fits for 
radiative cooling rates In a teaperature range of 0.002 KeV-100 



30 


KeV for 47 elenents 2 ^ Z £ 92. 

3 

Cooline rate i.e. cooling rate per free electron in erga/cm sec. 
is given by 

'■Z - i { 1®«10 [^. 

i=0 

where A(i) are available in [14]. 

3 

The energy loss rate in ergs/cn sec. can be found by 
multiplying the cooling rate by electron density and ion density. 

These fits include effects of all the three types of line 
radiations. 


2. 2. S Fusion Noubrons 

Fusion neutrons are the main source of loading among all 
sources considered, each time a deuterium and tritium fuses a 
neutron of high energy (14.5 HeV) is released. 

► ^He + ^n (14.5 MeV> 

These neutrons cause extensive damage to the first wall. 

Rate of this fusion reaction is 

R = njj . n^ < o . V >jj^ (2.44) 

From [16], the coefficient <o'.V>jj^ is given by 


<«'.V> 


DT 


3.68 X 10 
,2/3 


-12 


exp 


(KeV) 


f 19.94 

I ‘ :;t73 


CKeV) 


Jcm^/sec 


e e (2.45) 

So rate of neutron release from unit volume of plasma is given by 

(2.44), the total energy can be obtained by multiplying this 

number with energy of each neutrons (14.5 HeV). 
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2.3 Calcula'tlon algorithm tar local wall power loa<llng 

Becauae of the toroidal ayastetry and uniform porperty of 
plaama throughout the torua, (i.e. variationa of plaama 
parametera are aame in any part of the torua), it ia aufficient 
to conaider a amall thickneaa of the torua for thia analyaia. 
Reaulta obtained at any aection will be aimilar; Fig. (2.5>. 

Aaaume a amall thickneaa of unit length along toroidal 

direction ia cut out. Thia aection can be aaaumed to be a 
cylinder having elliptical croaa aection. Thia volume ia then 
divided into a number of cella, by drawing concentric ellipaea. 

The volume of each cell ia then calculated, uaing the 
following method. 

The area of the ahaded portion in Fig. (2.6) ia given by 
= tan ^ ^a/b ^ tand^ " j 


AA 


(2.46) 






33 


So cross sectional area of each of the cells can be determined hr 
subsequent subtraction from each other. Hultiplication of these 
surface areas with will qive volume of each cell. 

In section 2.2, analytical formulae are fiiven to compute the 
radiations in watt/cm^; when multiplied by volume of the cell, 
that will fiive total enerfiy radiated from that volume. 

Now, to compute the poloidal distribution of the radiation, 
which impinfies on the first wall the followinfi method is used. 

Select ine any arbitrary wall point and considerinfi small 
lenfith 'Al’ alonfi the poloidal direction of the wall in Fi£. 
C2.7), the fraction of the radiation eenerated in the cell that 
will reach A1 is 


„ A1 cos ft 

^ “ SiT osy 

where ft = ansle between line joinine the centre of A1 to 
of the cell and normal to Al. 


C2.47) 

centre 


ft - Anele CXBA) - Anfile (ZBC) C2.48) 

Ancle CZBC) = 90^, since CB is perpendicular to Al 


AmX. CXB») • Bl°P« CM) - alop. of (XTJ 

1 * slope of (AB) X slope (ZY) 

From each cell the total radiation can be calculated and 
usinc C2.48) and C2.49) the fraction of radiation that will reach 
the defined wall point can be calculated. 

So if '£’ s total radiation from any cell, 'A*. 

« fraction of 'E* which is reachinc 'B’ from 'A' * FxE . 
E^ = total radiation reachinc '&*,from all cells s s 




3S 


A 2 

So total enerey flux at 'B* = — ^ vatt/ca . 

A1 * A ^ 

2.3.1 Volumatrlc heat generation terns 

The local vail power loadings derived above is then used 
for tenperature solution of the first wall. Out of the various 
loading discussed, the breasstrahlung and neutrons causes 
voluaetric heat generation in the first wall. Breasstrahlung is 
high energy radiation and is not aono energetic, but distributed 
over a wide range of energies. The average photon energy coaing 
out as breasstrahlung radation is of the order of the electron 
teaperature [3]. This approxiaation is used here to calculate 
the aass absorption coefficient for the wall aaterial, 
corresponding to breasstrahlung radiation. 

The three aaln processes, which occur when gaaaa rays 
interact with aatter are- 

1. Photoelectric effect 

2. Coapton effect 

3. Pair production. 

Below 1.02 MeV of photon energy, the first two processes are only 
active [20]. 

Data for aass absorption coefficients for photon energy 

above 1 HeV are available in [21]. 

For photon energy below 1 HeV, the aass absorption 
coefficient can be calculated by approxiaate foraulas given in 
[ 20 ]. 


Photoelectric absorption coefficient. 
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t(E) ste 10“^^ N t } e“^*^ cm"^ (2.50) 

Compton absorption coefficient, 

<y(E) ste 1.25x10"^^ , ^ 5 ^ cm"^ (2.51) 

E |lofi^{2E/mc^>+ I j 

The high energy neutrons (~ 14.5 HeV) generated in the 

fusion reaction are sufficiently energetic to penetrate through 
the first wall. They cause most of the heating in the first 
wall, and damage the first wall. The neutron macroscopic 
cross-section (E) (which is equivalent to mass absorption 
coefficient in gamma rays) for different material are listed in 
[23]. 


The energy 

flux or intensity of 

beams 

(neutrons 

or fiamma 

2 

rays) I (watt /cm ) is equal to the product of 

flux 

and 

energy of 

the particles. 

For gamma rays I = 

^ (hv). 

hv 

s 

energy per 

photon, and in 

the case of neutrons. 

II 

AE, 

AE 

* energy 


loss per neutron reaction. 

In the case of monoenergetic photons or neutrons, I is 

directly proportional to so that I can be written as [24] 

I(x) = e"^* (2.52) 

I = intensity at x = 0. 
o 

Prom [24], the energy absorbed and consequently the heat 
generated in a particular location of an absorber where intensity 
of the radiation is I, is 

q"’ (X) = Ai I(x) (2.53) 

and specifically for neutrons 


q’'’(x) 


£ ^(x) AE 


(2.54) 
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where q"' = volumetric thermal aource strength (vatt/em )- 

The intensity variation with distance from (2.52) and 
(2.53) can be derived to be 

q*”(x) = q^' (2.55) 

where q^’ = volumetric thermal source strength at x=0. 

2.4 Tsmperature solution in the first wall 

In this section, two methods will be discussed for 
temperature solution in the first .wall, one for the surface 
heating (cyclotron, charge exchange neutrals, line radiation) and 
the other one for the volumetric heating terms (neutrons and 
bremsstrahlung) . 

Temperature solution for surface heating: 

Since a very small portion of the torus is considered for 
local wall power loading calculations, the first wall of that 
small portion can be assumed as a solid plate type element. A 
one dimensional heat flow condition is assumed. 

The heat conduction equation for this case (Fig. 2. 8) is the 
one dimensional polsson equation, 

2 ” ' 

^-4- + - — = ® (2.56) 

dx^ K 

K = thermal conductivity of the first wall material. 

Considering a thin layer of thickness Ax at distance x from 
midplane, the heat crossing plane x is given by 

- -KAg C2.S7) 

So 

q" 

- V * * 


t(x) 


S 


(2.58) 



Assuming the temperature at x » L (L = thickness) as 
t(x) = tj^ + ^ (L - X) (2.59) 

Temperature solution for volumetric heating: 

The case of flat slab is again treated. Fig. (2.9) shows a 
slab of thickness 'L' subjected to radiation from one side only, 
resulting in an exponential volumetric source strength 
distribution. Assuming one dimensional heat flow along positive 
X direction only, from (2.5£) 


= - K A 


using (2.54) 


dt 

Hx 


dt 

Zx 


e 


-jux 


(2.60) 


t(x) = - 

At X = L, t(x) = t. 




o 

e + c 


(2.61) 


which gives c = 


w 


* t. 


so t(x) = *1 + “ 15 ^ e - e J 


(2.62) 


Calculation of outer wall temperature: 

Coolant arrangement for the first wall is already discussed 
in section 1.5, i.e. coolant tube of the same material as the 
first wall are attached to the back of the wall and there is only 
one input and output of the tube. Hass flow rate and input 
temperature of coolant are known. 

Defining a coordinate 'a*,0 < a < 2na, on the perimeter of 
the coolant tube From general principle of convective heat 




Flc 2.8 Surface heating 



Fig 2.9 Volunetrlc heating 



40 


transfer , 


p C (T(s) - T. ) A u = 

P in xs 


heat load from plasma 
(2.43) 


between 0 to s 
where p - density of coolant 

Cp = specific heat of coolant 
A^^ = area of cross section of the coolant tube 

u = velocity of coolant. 

So T(e) = bulk temperature of the coolant at 's’ can be 
solved from (2.63). 

Consider inft a differential element of wall of unit area, 
h(Tj - T(a)) = plasma load at s (watt/cm^) (2.64) 
where h - heat transfer coefficient 

which can be obtained from the three dimensionless constamts. 
Reynold number, Nusslet number and Prandtl number 


Re » 

Dup 

P 

(Reynold number) 

(2.65) 





Pr = 

P 

K 

(Prandtl number) 

(2.66) 

Nu = 

0.023 

X X Pr”-* - 

(Nusselt Number) 



(2.67) 


where p = viscosity of liquid (coolant). 

Equation (2.64) fiives the outer wall temperature (t^) of the 
first wall. 

2.S Lind.'taHonB and Suggastlons 

1. It has been assumed that the average photon energy in 
bremsstrahlung radiation is of the order of electron temperature 



41 


of tho plasna. [3] . Since the breiBS8tr&hlun£ radiation ia one of 
the fltajor contributor to vail loading, a nore accurate method ia 
required to model the volumetric heatin£ cauaed by bremaatrahluns 
radiation by aplitting it in different energy £roupa. Hethod 
deacribed in Ref. [10] can be implemented for an accurate 
analyaia of the bremaatrahlunfi energy depoaition. 

2. For cyclotron emiaaion, it ia aaaumed that due to 
reflection from the wall and abaorption in the plaama, 98% of the 
cyclotron emiaaion ia loot, according to [3]. However thia ia 
only a rough eatimate. Tamor [3] did a tranaport analyaia of the 
balance between aynchroton-radiation and reabaorption for 
cylindrical and toroidal plaamaa and found that the plaama core 
radiatea atrongly, but the plaama ia actually heated near ita 
outer radiua, leading to a flatter temperature profile than would 
otherwlae exlat. Thua cyclotron radiation, tenda to cool the 
core of fuaion plaama and heat the outer periphery of the plaama. 
A full radiation tranaport calculation ia eaaential for an 
accurate eatimation of theae effecta. 

3. For the charge exchange neutrala aource term, it haa 
been aaaumed that number of charge exchange neutrala formed at a 
point, move tovarda the wall ia conatant throughout ita way upto 
the wall. Uhereaa due to further ionization and further charge 
exchange the number of neutrala formed at a point and the number 
of neutrala reaching the wall may differ. Detail calculation is 
required for thia attenuation of charge exchange neutrals. 

4. Only one dimensional heat conduction ia assumed in the 
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first wall without considering the toroidal shape of the first 
wall* i.e. by assuning the snail portion of the torus a cylinder 
and the wall is like a flat plate. Inclusion of toroidal 
geonetry (shape factor ) in the calculation of tenperature will 
give nore close value to the actual case. 



CHAPTER 3 

SYSTEM CODE FOR FIRST WALL TEMPERATURE CALCULATION 

3,1 A Brief Outline about the Code 

A code la developed for the first wall tenperature 
calculation of a tokanak reactor. The code is written in 
‘FORTRAN’ language. Hodel used for the purpose is discussed 
in Chapter 2. 

The code is written in a modular fashion, so that in 
future, as to requirement, with availability of better 
formulation, any part can easily be replaced. 

Some details about the code is given below: 

1. Title of the Code: First wall thermal analysis and 
temperature calculation in a tokamak reactor. 

2. Computer - H.P. 9000 
Installation place - l.I.T. Kanpur 

3. Operating system - Unix 

4. Language used - FORTRAN 

5. Number of lines — The main program is connected to 13 
link files (subroutines) and data files for polynomial fits. The 
entire code is *^850 lines (~ 30 Kb). 

6. Nature of physical problem - Calculation of all major 
types of radiations, with most of the parameters are inbuilt 
within the code, proper polynomial fits are used for this 
purpose, and based on that the one-dimensional temperature 
solution of the first wall, with a specific coolant arrangement 
system. 
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3 . 2 liqputf/'Dut.put of 'the Code 

Input to the 'FORTRAN' code includes the followlns: 

The corresponding variables used in the code are shown in 
bold capitals. 

1. Plasna dimensions 

a) Plasma major radius (in cm) *XNJR* 

b) Plasma minor radius (in cm) *XIINR* 

2. First wall dimensions 

a) Major radius of first wall (in cm) ‘XlOiaf* 

b> Minor radius of first wall (in cm) *XilMRW* 

3. Density and temperature profiles 

The code will assume a flat type density and temperature 
profile, if *XC(pnON* La set to any number less than 2. For any 

other value of *XOPTION* , the code will go for a peaked profile. 

The constant oi , a. of profile type, eqn. (2. 1,2. 3) can be varied, 

XI X 

by changing 'ALPHAN* and ‘ALPHAT* . 

4. Impurity contents 

a) Ratio of, density of hydrogen to density of electron 

•HYDIMP* 

b) Ratio of, density of carbon to density of electron 

'CARBIMP* 

5. Magnetic field strength 

a) Magnetic field strength in 'Gauss' ‘BMAe* 

6. First wall material 

a) Thickness of first wall material (cm) ‘WALLTHICK* 

b) Conductivity of first wall material (W/cm C) •THCWALL* 
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c) Neutron removal croee-eection for the first wall 
material *XMUE C^ble variable appears in Link~‘12). 

7. Coolant tube 

a) Thickness of coolant tube wall (cm) •PIPETHICaC* 

8. Coolant 

The following input informations are required for the 
coolant. They are dependent on coolant system pressure. 

a) Density of coolant (kg/cm^) *DG14C(XX<* . 

b) Viscosity of coolant (kg/cm-sec) ‘VISCOSITT* 

c) Thermal conductivity of coolant (U/cm ®C) •TIKMOOL* 

d) Inlet temperature of coolant C®C) *TSIIPINCOOL> 

9. Uall point position, where the temperature profile is 
required in degree. *THETAW* , THETAU = 0 indicates minor axis. 

The code performs the following calculations: 

1. The spatial profiles of the intensities of line and 
recombination radiation, charge exchange, neutrals, 
bremsstrahlung, cyclotron radiation and fusion neutrons. 

2. Fraction of particle flux that hits the wall point of 
interest . 

3. Based on item 1, calculation of the surface and 

volumetric heat load on the wall as a function of poloidal angle. 

4. Calculate the heat transfer coefficient and hence the 

temperature profile through the wall. The temperature profile at 
50 points of the wall are redirected to an output file 


•TEMP.DAT’ . 
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3. 3 Flow Chart 
















CHAPTER 4. 
RESULTS 


Ttid codo WAS ‘tos'tsd with ssvsrAl ssts of inpuf p&rAnstsrs* 
includinfi the 'ADITYA’ parameters. The results shown here are 
for the followinfi 'ADITYA' parameters. Other parameters, whenever 
used is indicated separately. 

1. Electron temperature = 15 KeV (max) 

2 . Ion temperature = 15 KeV (max) 

3. Electron density = 4el4 f cm“^ (max) 

4. Ion density = 4el4 # cm (max) 

5. Hafinetic field strenfith « 1.5 T 

6. Plasma minor radius s 30 cm 

7. Plasma major radius « 60 cm 

8. First wall minor radius » 32 cm 

9. First wall major radius = 62 cm 

10. Carbon density in plasma = 0.1 x electron density 

11. Bydroflen density in plasma = 0.1 x electron density 

12. Thickness of first wall = 2 cm 

13. First wall material = S.S. 314 

14. Coolant tube material ~ S.S. 314 

15. Coolant tube wall thickness » 0.5 cm 

16. Coolant tube diameter = 1.0 cm. 

17. Hass flow rate of coolant = 90.7E-04 Kg/cm^- Sec 

18. Coolant is ordinary water. 

19. Inlet temperature of coolant = 20®C 

The poloidal distribution of the total flux at the inside 
surface of the first wall (facing the plasma) both for flat 



profile and peaked profile is shown in fig 4.1 and fig. 4.2 
respoctively.lt is seen from the curve that the loading is 
maxinuro at e - 90 poloidal angle) and 6 = 270 for both 

type of profiles. The distribution of the flux is ssrnnetrical 
about the major axis. 

Flat type profile: 

2 

naximum value of flux » 86.150 watt /cm at d = 90 & 270 
Hinimum value of flux = 49.498 vatt/cm^at d » 0 & 180 
Peaked type profile : 

2 

Maximum value of flux = 67.623 vatt/cm at 6 =90 & 270 

2 

Minimum value of flux = 34.277 watt/cm at d = 0 & 180 
The poloidal distribution of temperature at the in side 
surface of first wall ( facing the plasma ) for both flat and 
peaked profile is shown in fig 4.3 and fig 4.4 respectively. As 
expected from fig 4.1 and fig 4.2 the maximum temperature is at ^ 
= 90 and 270 and the minimum temperature is at d = 0 and 180 . 
Flat typ* profila i 

Maximum temperature = 832.682®C at d = 90 a 270 

Minimum temperature = 505.008^C at d = 0 & 180 

Paakad typa profila i 

Maximum temperature = 670.600^C at d = 90 a 270 

Minimum temperature = 371.523®C at d = 0 a 180 

The temperature profile of the first wall at one of the 
hottest point ( d = 90 ) is shown for the flat and peaked type 
profile in fig 4.5 and fig 4.6 respectively. 


Flatr type profile i 

At ( X = 0 cok ) aurface facing the plaaaa Tenp.= 832. f 82 

At ( X = 2 cm ) surface away from the plasma Temp.s 568.366 

peaked type profile i 

AT ( X = 0 cm ) surface facing the plasma Temp. = 670.600 

At ( X = 2 cm ) surface away from the plasma Temp. = 459.836 

Fig 4.7 to fig 4.10 shows the variation of source term 
strength with electron temperature. Fig 4.11 shows the variation 
of the source term (charge exchange neutrals) with ion 
temperature .These results are for the same parameter set given 
above except the Impurity levels. Uhlch Is assumed as : 

The ratio of density of carbon to density of electron In the 
plasma 0.01 

The ratio of density of hydrogen to density of electron In 
the plasma « 0.01 

Fig 4.7 shows the variation of bremsstrahlung radiation (watt/cs^ 
) with the electron temperature. As expected from eq.2.18 the 
variation Is propertlonal to square root of electron temperature. 
Fig 4.8 shows the variation of cyclotron radiation (watt/ci? ) 
against the electron temperature. As expected from eq.2.33 the 
variation Is linear. 

Fig 4.9 shows the variation of line radiation with electron 
temperature. Line radiation la low as expected, because at these 
high temperatures Impurities are almost fully Ionized, l.e there 
e.re no bound electrons. In reality there will be a small radiating 
region near the edge where the temperature Is low. 




POLOIDAL ANGLE 


Fig. 4.1 Poloidal diatrlbtttlon of total flux at th« inside 
surface of first wall (flat profile) 






TEMPERATURE (C) 



POLOIDAL ANGLE 


Fig. 4.3 Poloidal distribution of oaxlaun teapsrsturs 
inside surface of first wall (flat profile) 
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Fig 4.10 shows th® varlstlofi o£ ■ nsutron radlatloa with ths 
electron temperature. It Is seen from the figure that neutron 
source term Increases with Increase of electron temperature. That 
Is because the rate of fusion reaction Increases with 

temperature. Also this Is seen that among all the source terms the 
maximum contribution Is coming from the neutron term. 

Fig 4.11 shows the variation of charge exchange neutrals source 
term with the Ion temperature. This source term Is very low as 
expected , because at such high temperature the recombination 
reaction rates which are responsible for neutrals formation is 
small. Since It Is assumed that energy carried by the neutrals 
Is 1.5 times the ion temperature the curve shows a increase 
tendency with increase In ion temperature. 

The code Is tested with variation of some important 

3 

parameters and corresponding source terms (watt/cm ) Is shown in 
table 4.1 to table 4. 4.. Only the shown parameter Is varied and 
other parameter for these results are as given for "ADITTA”. 

Table 4.1 shows the density variation of electron (Ion density Is 
same as electron density ) .corresponding strength of source 
terms are tabulated. As expected all the terms show an Increase 
In value with increase In density. 

Table 4.2 shows the variation of electron temperature and 
corresponding strength of source terms .Ulth Increase In electron 
temperature all the source terms strength Increases. Table 4.3 
shows the variation of Impurities level. Except cyclotron all 
other terms are affected. Since with Increase ’of Impurity the 
amount of D an T reduced l.e. because electron density and Ion 



62 


density is fixed, so neutron contribution which is dependent on 
the density of D and T reduces. 

Table 4.4 shows the variation of nagnetic field strength and it 
is seen that only cyclotron radiation is affected by this which 
was expected. 
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TABLE 4.1 

VARIATION OF THE SOURCE TERNS UITH THE ELECTRON DENSITY 


E DENSITY 

BR 

CY 

CZ 

LINE 

NT 

0.20E+14 

0.107E-02 

0.838E-03 

0.268E-08 

0.375E-03 

0.373E-01 

0.40E-I-14 

0.430E-02 

0.167E-02 

0.107E-07 

0.150E-02 

0.149E+00 

0.60E-I-14 

0.968E-02 

0.251E-02 

0.241E-07 

0.337E-02 

0.336E+00 

0.80E-I-14 

0.172E-01 

0.335E-02 

0.430E-07 

0.600E-02 

0.598E-t-00 

O.lOE+15 

0.269E-01 

0.419E“02 

0.672B-07 

0.937E-02 

0.934E-)-00 

0.12E+15 

0.387E-01 

0.503E-02 

0.967E-07 

0.135E-01 

0.134E-i^01 

0.14E+15 

0.527E-01 

0.586E-02 

0.131E-06 

0.183E-01 

0.183E+01 

0.16E+15 

0.688E-01 

0.670E-02 

0.172E-06 

0.240B-01 

0.239E-i-01 

0.18E+15 

0.871E-01 

0.754E-02 

0.217E-06 

0.303E-01 

0.302E-i-01 

0.20E-I-15 

0.107E+00 

0.838E-02 

0.268E-06 

0.375E-01 

0.373E+01 


E.DENSITY»Electron density (f . /cs^) 

3 

BRs Brensatrahlung radiation (Uatt/cs ) 

3 

C7s Cyclotron radiation (Uatt/cn ) 

3 

CX= Charge Exchange neutrals (Uatt/cn ) 
LINES Line radiation (Uatt/cn^) 

3 

NTs Neutron radiation (Uatt/cn ) 
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TABLE 4.2 

VARIATION OF THE SOURCE TERNS UITH ELECTRON TENPERATURE 


E.TEMP, 

BR 

CY 

CX 

LINE 

NT 

0.20B+04 

0.157E+00 

0.223E>02 

0.444E-06 

0.898E-01 

0.249E-01 

0.40E+04 

0.222E+00 

0.447E-02 

0.536E-06 

0.974E-01 

O.dlOE+OO 

0.60E-i-04 

0.272E+00 

0.670E-02 

0.614E>06 

O.lOBE-i-OO 

0.345E+01 

0.80E+04 

0.314E+00 

0.894E-02 

0.692E-06 

0.118E-I-00 

0.345E-I-01 

O.lOE+05 

0.351E-I-00 

O.lllE-01 

0.779E-06 

O.lXBE-i-OO 

0.608E+81 

0.12E’i-05 

0.385E+00 

0.134E-01 

0.881E-06 

0.137E+00 

0.929E-f-01 

0.14E+05 

0.41SE+00 

0.156E-01 

O.lOOE-05 

0.146E+00 

0.129E+02 

0.16E+05 

0.444E-I-00 

0.178E-01 

0.115E-05 

0.153E-)-00 

0.170E+02 

0.18E+05 

0.471E+05 

0.201E-01 

0.133E-05 

0.161E+00 

0.213E-i-02 

0.20E+05 

0.497E+00 

0.223E-01 

0.156E-05 

0.168E+00 

0.258E+02 


E.TBHP. “Electron temperature (ev) 

3 

BR“Br ernes trahlung radiation (Uatt/cm ) 

3 

cy=Cyclotron radiation (Uatt/cm ) 

CX“ Charge Exchange neutrals (Uatt/cm^) 
LlHE=Litte radiation (Uatt/cm^) 

3 

NT* Neutron radiation (Watt/cm ) 
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TABLE-4. 3 

VARIATION OF THE SOURCE TERMS HITH MAGNETIC FIELD STRENGTH 


BMAG 

BR 

CY 

CX 

LINE 

NT 

O.lOEtOB 

0.430E+00 

0.745E-02 

0.107E-05 

0.150E+00 

0.199E+02 

0.20E+05 

0.430E+00 

0.298E-01 

0.107E-05 

0.i50E+00 

0.149E+02 

0.30E+05 

0.430E+00 

0.670E-01 

0.107E-05 

0.150E+00 

0.149E+02 

0.40E’t-0S 

0.430E+00 

0.119Et00 

0.107E-05 

O.lSOE-i-OO 

0.149E+02 

0.50E+05 

0.430E-I-00 

O.lSfE+OO 

0.107E-05 

0.150E+00 

0.149E+02 


BHAGsHaenetlc field atrength 

3 

BRsBrenaatrahlung radiation (Uatt/cn ) 

3 

CY“Cyclotron radiation (I7att/ca ) 

3 

CX=Charge Exchange neutrals (llatt/cn ) 
LINE=Llne radiation (Uatt/cn^) 

3 

NT^Neutron radiation (llatt/ca ) 



66 


TABLE A. 4 



VARIATION 

OF SOURCE 

TERNS WITH 

IMPURITY 

CONTENTS 


HIMP 

GIMP 

BR 

GY 

CX 

LINE 

NT 

0*1E-01 

O.lE-03 

0.33E+00 

O.KE-Ol 

0.103E-05 

0.15E-02 

0.16E+02 

O.lE-02 

O.lE-02 

0.34E+00 

0.14E-01 

0.107E-07 

0.17E+02 * 

0.17E+02 

O.lE-03 

O.lE-01 

0.43E-t-00 

0.16E-01 

0.107E-09 

O.lSE-i-OO 

0.15E+02 


HinP^Hydrogen denaity /Electron density 
ClNP»Carbon density /Electron density 

3 

BR^Brensstrahlung radiation (Uatt/cm ) 
CYaCyclotron radiation (Watt/ca^) 

3 

CX-Charee Exchanged neutrals (tfatt/ca ) 
LINE=Line radiation (Uatt/cm^) 

3 

NT»Neutron radiation (Uatt/cm ) 
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